Ion- Wake induced anomaly of Dust Lattice Mode in presence of external 

Magnetic field 

Saurav Bhattacharjee la \ Nilakshi Das lb ) 
Department of Physics, Tezpur University, Tezpur, Assam, 784 028, India 

o" 

Email: sauravtsk.bhattacharjee@gmail.com 

< 

(N 

we report a theoretical investigation of DL mode (Dust Lattice) in 2D Yukawa 

&\ 
♦pi! crystal in presence of asymmetric ion flow and external magnetic field perpen- 

g: 

c« ' dicular to the crystal plane. Two mutually perpendicular modes are found to 

'Si- 

_o 

cal direction is strongly affected due to the formation of ion-wake. This causes 
Q^ anisotropy in interaction strength along two mutually perpendicular directions. 

Both hybrid modes are studied as characteristics of different ion flow speed and 
CN ' magnetic field strength. The study shows a fluctuation in DL mode frequency 
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be coupled due to Lorentz force. Interaction among the dust grains along verti- 



driven by the strength of the particle-wake interaction. Effect of ion flow on po- 
larization of the hybrid wave amplitudes is discussed in details. Results show a 
possible mechanism of anomalous phase transition in dusty plasma. 
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1. Introduction 

The presence of micron sized solid particles in plasma adds a new dimension to it. Be- 
cause of large and fluctuating charge on dust grains various novel and rich phenomena arise 
in plasma. Highly charged dust grains of nanometer to micrometer size interact with each 
other through repulsive Yukawa potential. Under the influence of strong repulsive poten- 
tial, particles may form fluid and crystalline structures. This provides an interesting phase 
of dusty plasma that facilitates the study of various problems of fundamental physics like 
phase transition, crystal defects etc. Research in dusty plasma has observed an exponential 
growth after the discovery of plasma crystal by two pioneering groups led by Chu et. al. [1] 
and Thomas et.al.[2] in 1994. The crystal structures include, starting from the ID chain, to 
2D and 3D cubic crystals having different lattice geometries as Triangular, FCC, BCC and 
Hexagonal etc. Recent experiments in microgravity conditions provide an extended system 
of complex plasma for further investigations. 

The discovery of low frequency dust mode in plasma [3] has triggered interest in study- 
ing other waves like dust ion-acoustic wave(DIA)[4] , electrostatic dust ion cyclotron waves 
(EDIC)[5], dust lattice waves (DLW) [6] . Melands0[7] predicted theoretically the existence 
of dust lattice waves in 1996. This mode was experimentally observed in low temperature 
laboratory dusty plasma discharges by Homann eta/. [8] in 1997. In recent years, a variety 
of linear and nonlinear modes of the dust lattice waves have been investigated both exper- 
imentally and theoretically[9-10] . In presence of magnetic field perpendicular to crystal 
plane, two mutually perpendicular modes of a 2D yukawa crystal found to be coupled due to 



lorentz force [11]. This leads to two type of hybrid modes of high and low frequency branches 
instead of two independent modes. 

As dust crystals are formed near the sheath region, the plasma flow dynamics near the 
sheath give rise to various attractive potentials such as shadowing and wake potential etc. in 
addition to the repulsive Yukawa type of interaction among dust grains. Closer to the lower 
electrode the velocity of the ion stream towards the negatively charged electrode may lead to 
instabilities like fluctuation in particle positions [12- 13] and develop a new kind of potential 
named as wake potential. This instability leaves an impact on the dust lattice modes as 
well. It is identified as the possible cause for melting of dusty plasma crystal. Ivlcv ei.aZ.[14] 
have given an interesting theory on anisotropic dust lattice mode. They have also reported 
the decrease in frequency of the dust lattice modes due to particle wake interaction. Joyce 
e£.a/.[15] in 2001 for the first time presented a simulation model to study the behaviour of 
plasma crystal, incorporating both short range Yukawa potential as well as long range wake 
potential. They reported the formation of vertical strings of dust grains along the downward 
flow direction. 

Ion- wake structure and their dynamics have been extensively studied in absence of mag- 
netic field. In presence of magnetic field, very limited theoretical and experimental studies 
[16-17] are available. In the limit of linear response theory, Nambu et.al. [16] have reported 
damping of wake strength in presence of magnetic field. This leads to a weak effective inter- 
action strength between dust grains along vertical direction as compared to unmagnetized 
case and is verified experiment ally [17] in recent year. In a recent paper Saurav et.al. [18] 



have investigated the effect of wake potential on Coulomb crystallization in presence of mag- 
netic field. They observed a reduction in RDF (Radial Distribution Function) peak heights 
with increased strength of wake potential and reported a ion-wake induced dislocation of 
dust grains from regular crystal arrangement. This modification in particle- wake interaction 
in presence of magnetic field with the strength of the ion flow and its effect on various eigen 
modes is an important issue in complex plasma. 

The phase behaviour of fluids depends solely on the intermolecular interaction. Un- 
derstanding this mechanism is one of the biggest challenges in statistical thermodynamics. 
Ion-wake formation and its dynamics found to play a vital role in understanding the forma- 
tion of coulomb crystals and their phase transition in magnetized complex plasma. In this 
letter we presented a theoretical study to investigate the behaviour of Dust lattice modes in 
presence of magnetic field perpendicular to crystal plane and ion flow. Polarization of both 
the hybrid wave amplitudes and their interplay with the magnetic field have been discussed 
in details. Results show a possible kind of anomalous phase transition [19] in dusty plasma. 

2. Theoretical Formulation 

2.1. Interaction Potential among dust grains 

We consider a 2D uniform dusty plasma with stationary heavier dust grains in presence 
of an external magntetic field B D x perpendicular to the downstream ion flow. Under the 
influence of magnetic field, the ion flow dynamics get affected and subsequently the dielectric 
response function [16] of the plasma medium also changes. In presence of low frequency 
wave (w-CWce) ,with phase velocity smaller than electron thermal velocity, electrons will have 
a Boltzmanian distribution. In normalized scale of electron Debye length Xoe, the dispersion 
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relation can be written as: 
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k 2 e(u, k) " k 2 + 1 (jfc2 + l)M 2 (fc|| 2 - A; R 2 )(A;|| 2 - k Im 2 ) 
where it is assumed that k\y^>k y and k 2 =k z 2 +k y 2 , k z (\\) and k y (±.) are wave vectors 
parallel and perpendicular to the direction of ion flow. Here M=^ , is the mach number 
having Cj=Ai) e u;™ as ion acoustic wave velocity and Uj is the streaming velocity. / = -^ 
Upi and u C i are the ion plasma frequency and ion cyclotron frequency respectively. k# and 
kj m corresponds to real and imaginary screening modes as a function of kj_, in a complex kn 
plane. 



W = -JkJ + 1 - _ - _) ± ^( ti > + !-_-_)* + ^ (2) 
In the limit of kj_>jjj2 , magnetic field B<1.0T, for all M, both the pole can be approximated 



as k R ttJj±^, k lm xi±iyjk± 2 + 1 - jjja-. 

Using the dispersion relation in equation (1) one can write the interaction potential 
between the test dust particulates in presence of magnetic field as 

»<»•" = (2^/^w'ftt + & + ijiiw -lw -o ' e " J (3) 

where v=\Jy 2 + z 2 . Integrating the first part on the R.H.S of the above equation gives 

Yukawa potential 

$ D = ^-exp(^) (4) 
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Integrating the second part on R.H.S of equation (3) for the real pole gives oscillatory wake 



potential 
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In the limit of / 3>1, above equation shows a damping in the strength of the wake potential. 

2.2. Dispersion relation of Dust Lattice Wave (DLW) 

We have considered a simplified model of monolayer hexagonal lattice in Y-Z plane with 
ion flow along downward z-direction and magnetic field perpendicular to the crystal plane. 
The lattice mode propagate along Y-direction parallel to the fundamental lattice vector 
directed from center of the cell to one of its nearest neighbors. The equation of motion along 
two transverse directions are given by [11] 

d 2 y n i v^ r , / v , „ „ dz n 
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Yl K vn m (ym - Vmi) + QdB —^- (6) 
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d 2 z nl ^ ,, %, 



XI [■ K *»m + ^«o„J(^l-«ml)-Q < i-B -^ (7) 



t,2 / , , L"2nm ' "ZiUnmJ \~«± ~mi; -^a~o ,, 

m=l,n7^m 



Where yi n ,m)i an d z {n,m)i are the displacement of the particles from their equilibrium posi- 
tions. K( y ,z) nm are the spring constants corresponding to Yukawa type of interaction along 
y and z direction and K ZWnm is the spring constant corresponding to wake potential along 
vertical z direction. 

Using equation (4) we can write the normalized spring constants K y , z corresponding to 
Yukawa potential along y and z direction as 

d 2 $> D Texp(-V2K) r- 2 

Ky > z = d^z? ly ' Z=a = 2W l ] {) 

Here a is average interparticle distance, T is coulomb coupling parameter and k = y— stands 

for screening constant. 
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Similarly from equation (5)we get the spring constant for the wake potential as 
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zw M(} + M 2 )(} + M 2 -l) {) 

This leads to the effective spring constant along z-direction as 

K z = K z (Yukawa) + K zw (Wake) (10) 

As we have considered that DL mode propagates along Y direction, the oscillations of the 
wave modes are y n i=A y exp[i(k y a — cot)] and z n i=A z exp[i(k y a — cut)] respectively along two 
mutually perpendicular directions. Substituting y n \ and z n \ in equation (6) and (7) we get 
the dispersion relation as: 

N N 

[u 2 + Yl K y nm ^M ik v a ] ~ 1 )}l uj2 + Y ( K *nm + K ZWm )(exp[ik y a] - 1)] - w c V = 

(11) 

The dispersion relation indicates that two transverse modes get coupled by Lorentz force 

and there exist two hybrid modes instead of independent modes. 

In the limit of k y a<^l, cos(k y a) ~ 1 — j(k y a) 2 + .... equation (11) gives two roots as 






:+ Yl m ynm +K Znm + K ZWm )sm 2 {^-), [Upper Hybrid] (12) 
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Co> 



(13) 



Here the frequency is normalized by ujn = \ B vi d , where Kb, Tj, ma and A^ be the 

y m d A o e 

Boltzmann constant, dust temperature, dust mass and dust Debye length respectively. 



Substituting the displacements y nm and z nm in equation (6) and (7) we can get the wave 
amplitude ratio along y and z direction as: 



Unm . I,, u 2 + 2(K z + K zw )(cos(k y a) - 1) 

- ifi, where \fi\--=\ " (14) 

ur + 2Ky(cos(kya) — 1) 



Using equation (14) one can write the polarization P of the wave amplitudes as 
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(15) 



P=l corresponds to pure longitudinal mode and P=0 corresponds to pure transverse mode. 
3. Results and Discussion 

The Debye cloud surrounding the dust particles gets distorted by the downward ion flow 
and becomes elongated in lower vertical (z-direction) direction. Amount of distortion and 
the scale length of elongation depends mainly on the strength of the magnetic field and flow 
velocity of the downward ion stream towards the plasma sheath. As there is no attractive 
type of potential in the horizontal direction, this excess positive charge on the lower side of 
the Debye cloud attracts the dust particles just below it to align along vertical direction. This 
process of alignment has a vital dependency on the ion-wake induced vertical confinement 
frequency. It is obvious from equation (5) that frequency of the wake potential increases 
with decrease in Mach number for a given magnetic field strength. As dust particles are 
confined at the trail of the ion-wake, alignment continues up to a certain value of confine- 
ment frequency and beyond that dust particles are no more confined and they restore their 
original configuration once again, as demonstrated in Fig. 1(a) . Fig. 1(b) clearly shows the 
transition of vertical confinement strength with ion flow speed. 



FIG. 1: (color online) (a) Vertical confinement mechanism of particle- wake interaction with 
increase in strength of wake potential from left to right, (b) Effective spring constant K z vs. 
M(mach number). 

Due to this interaction between dust particles and ion-wake, partciles are first dislocated 
from their regular orientation and then again restore their original configuration. This entire 
process leaves a strong impact on the dust lattice modes. It is visible from the dispersion 
relations (12) and (13), shown in Fig. 2, for different magnetic field strength and ion flow 
velocity. 

FIG. 2: (color online) The plots show the dispersion relation profiles with different Mach 
number M and magnetic field strength B for coulomb coupling parameter T=824 and screen- 
ing strength k—2.9. (a) Shows the normalized frequency dependence of both Upper Hybrid 
(UH) (solid line) and Lower Hybrid (dotted line) waves with M=0.3(green), 0.18(blue) and 
O.ll(red), for a given magnetic field B=0.08T. (b) M=0.6 (green), 0.4(blue) and 0.25(red) 
with B=0.06T. (c) M=0.5 (green), 0.3 (blue) and 0.15(red) with B=0.04T. 

It is evident from equations (14) and (15) that polarization of the hybrid wave amplitudes 
varies with ion flow speed as well as the magnetic field strength. At comparatively higher 
ion flow speed (M>0.4), the strength of the wake potential is not sufficient enough to affect 
the polarization of the wave amplitudes as can be seen from FIG 3(a). Above discussion 
clears the fact that effective repulsive strength along vertical direction first decreases with 
decrease in M value (strong confinement zone) due to enhancement in strength of the wake 
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potential. In that region we find the lower hybrid mode becomes more longitudinal in nature 
as compared to upper hybrid for a given magnetic field strength and wave number. 

FIG. 3: (color online) The profiles show (a) Polarization P vs. Mach number M, for 
both upper hybrid (UH) (green solid) and lower hybrid (LH) (green dotted) waves for a 
given wave number ka/7r=0.6 and magnetic field strength B=0.08T. (b) P vs. M for a given 
wave number ka/7r=0.6 and different magnetic fields B=0.08T, UH (green solid), LH (green 
dotted). B=0.06T, UH (blue solid), LH (blue dotted), (c) Strength of spring constant k 2 vs. 
f(B)=uj p i 2 /u C i 2 for a given plasma frequency 0^=1.86 xlO 7 s _1 and Mach Number M=0.5. 

On further decrease in Mach number (weak confinement zone) below M< 0.24, upper 
hybrid becomes more longitudinal in nature as compared to lower hybrid, shown in Fig.3(a). 
It can be seen from Fig. 3(c) that effective spring constants along vertical direction de- 
creases with increase in strength of the magnetic field (decrease in / value) and agrees with 
the experimental result [17]. This causes shifting of the critical Mach number, from which 
confinement starts dominating, to a lower value, with decrease in strength of the magnetic 
field as shown in FIG. 3(b). This ion- wake induced periodic transition in crystal phase from 
a ordered to disorderd state and again get back to the ordered state showing an analogy 
with the anomalous phase transitions in complex fluids. 

4. Conclusions 

In the present piece of work we have elaborately discussed the effect of wake potential 
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on dust lattice modes in terms of ion flow speed and external magnetic field for a 2D crystal 
structure. The observations can be summarized as follows: 

(a) The ion flow near the sheath region gives rise to attractive wake potential along 
vertical direction and this causes anisotropy in the interaction strength along two mutually 
perpendicular directions. 

(b) It is demonstrated in FIG. 1 that with increase in strength of the wake potential, 
the dust particles start aligning along vertical direction from their regular orientations in 
the crystal plane. On further increase in strength, confinement stops and particles regain 
their original configuration. Impact of this mechanism on DL modes is clear from dispersion 
relations in Fig.2. 

(c) Polarization of both the hybrid wave amplitudes is also affected by the strength of 
the wake potential. Within the strong confinement zone upper hybrid mode is found to be 
more transversely polarized as compared to lower hybrid for a given magnetic field and wave 
number. On further decrease in M value, in the weak confinement zone of interaction, po- 
larization of both the hybrid wave amplitudes reverses in character and upper hybrid mode 
becomes more longitudinally polarized as compared to lower hybrid. It can be seen from 
FIG. 3(a). 

(d) With increase in strength of the magnetic field, effective interaction along z-increases. 
As a result, the critical Mach number is shifted to a lower value with decrease in strength 
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of the magnetic field for a given wave number. It is shown in FIG. 3(b). 

(e) Finally, we can say the particle-wake interaction is responsible for decrease in fre- 
quencies of the DL modes up to certain strength of wake potential. Beyond that DL mode 
frequency increases again with the strength of the wake potential. 
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